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Making a soft relativistic mean-field eqnation of state stiffer at high density 
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We study relativistic mean-field (RMF) models including nucleons interacting with scalar, vector 
and iso-vector mean fields and mean-field self- and cross-interaction terms. Usually, in such models 
the magnitude of the scalar held increases monotonically with the nucleon density, and the nucleon 
effective mass decreases. We demonstrate that the latter quantity stops decreasing and the equation 
of state stiffens, provided the mean-held self-interaction potential rises sharply in a narrow vicinity 
of the values of mean helds corresponding to nucleon densities n n« > no, where no is the nuclear 
saturation density. As a result the limiting neutron star mass increases. This procedure ohers a 
simple way to stiffen the equation of state at densities above n* without altering it at densities 
n ^ no- The developed scheme allows a neutron star application of the RMF models, which are 
well htted to hnite nuclei but do not fulhll the experimental constraint on the limiting neutron star 
mass. The exemplary application of the method to the well-known FSUGold model allows us to 
increase the limiting neutron star mass from 1.72 Mq to M > 2.01 Mq. 

PACS numbers: 21.65.Cd, 26.60.-c, 


A relativistic mean-field (RMF) model proposed and 
advertised in [H-Q is a convenient vehicle for construction 
of the equation of state (EoS) of baryon matter, which 
preserves causality. Various RMF models are success¬ 
fully used to describe neutron stars (NS) and heavy-ion 
collisions, see [1-11]. The models prove to be also ap¬ 
plicable to atomic nuclei, where a high precision in the 
description of gross properties of finite nuclei close to the 
valley of stability can be achieved 0, [Ij . An accurately 
calibrateclparametrization of the RMF model was intro¬ 
duced in 01, known as the FSUgold model, which allows 
one to compute the ground state properties of finite nu¬ 
clei and neutron rich matter. 

In order to describe finite nuclei well, one includes 
the self- and cross-interaction terms of meson fields 0 - 
ini and/or nonlinear derivative couplings (ll |. However 
these models yield a rather soft EoS and cannot describe 
the heavy NSs with masses of (2.01 ± O.O 4 )M 0 , where 
Mq is the mass of the Sun, which were recently unam¬ 
biguously identified experimentally Bin. Additional 
complication arises if hyperons are included in the consid¬ 
eration, since in their presence the EoS softens even more, 
cf. [101 ■ To reconcile the appropriate description of the 
properties of atomic nuclei and the NS mass constraints, 
one exploits density-dependent coupling constants 0 - 
Similarly, one could use mean-field-dependent coupling 
constants, cf. [110. 

In this Rapid Communication, we will demonstrate 
that, if the mean-held self-interaction potential rises 
sharply in a narrow vicinity of the values of mean helds 
corresponding to nucleon densities n > > no, where 

no is the nuclear saturation density, the nucleon ef¬ 
fective mass saturates and the EoS stiffens. As a re¬ 
sult the limiting NS mass may increase above the value 
(2.01 ± 0.04)Mq. This procedure offers a simple way to 
stiffen the EoS at high densities without altering it at 


densities n < no. 

For illustration purposes we work with the standard 
non-linear Walecka (NEW) model 0 dehned by the La- 
grangian 


C = (idu - - 9ptPu)u^ - mN + gacr n 

topii, — dfiiOu djytOfj, , PfiL' — ^pPu ^uPp, ■ (1) 


Here 'Fa? = is the isospin doublet of Dirac 

bispinors for protons (p) and neutrons (n) with the bare 
mass rriN = 938 MeV; a, w, and p denote the helds of 
scalar, vector, and isovector-vector mesons with masses 
ni(T, TOtj, and nip, respectively; 7 ^ {p — 0,1, 2 ,3) are the 
standard Dirac 7 matrices; and t stands for the operator 
of the baryon isospin. The £int may include the scalar 
and vector-held self-interaction terms and the held cross¬ 
interaction terms. Following [I0| we will hrst use the 
simplest form of the £int reducing to the scalar-held self¬ 
interaction potential Ant = —U{a) = —buiN {ga-o')^/0 — 
c (<?crCT)‘^/4; g^, ptj, gp, b, and c being coupling constants. 

Equations of motion for the mesonic helds are treated 
in the mean-held approximation with the solutions 


a XX dpt't / \ X 9wN / N 

Pp — baSOpQ— 2 ^n), tOp — OpO ^ (^p “f tin) , 

Zrrip 

mla+ U'{a) = ga{ns,p + ns,n), (2) 


where ns,i, i = n,p, is the nucleon scalar density, which 
can be expressed through the nucleon density rii as 

ns,i = 0n^[xi{xf + 1)^/^ - log((x/ -I- 1)^^^ -I- Xi)]/(2xi), 

where Xi = p-p,i/ra*jq, and p-pp = ( 37 r^ni)^'^^ is the Fermi 
momentum. The effective nucleon mass = mAr(l —/) 
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depends on the a field expressed through the dimension¬ 
less variable / = g^jO jraj^. Then the energy density fol¬ 
lowing from the Lagrangian m is given by 


Ei^Tlp., Tln^ /*) — -£^kin(^p; '^ni “t“ Eyijlp-, ‘^n) E , 

PFn PFp 

^^kin(np,n„,/) = (j + j + 

0 0 

Ey ijlp 1 ) ~Z ^ ij^p E Tiji ) “h — ^ {P'p ) ? 

ZTTIj^ oTTijy 


EAf) = 


m%P 

2Ci 


+ U{f), m*j^{f) = mN{l-f), (3) 


Cj = gjiriN/rrij with j = cr, w, p. In terms of the / 
variable U{f) = m%{bf^/3 + cf^jA). The pressure is 
given by: 


P{j^n: ^p: f') — I^kin (^n : f')A~ Py ij^n: ^p') T Pa (^) ; 

PPn PPp 


fkin(f 


n ? '^p 


) /) — o( / + / ) 


p^dp 


P 


VP' 




Pvinn.rip) = Evinp,n„), P„{f) =-Ec{f). (4) 


Five model parameters Cuj, Ca, Cp, b and c are tuned 
to fit the values of the saturation density ng, the bind¬ 
ing energy per nucleon Eg = £'(no/ 2 , no/ 2 , /o)/no - tun, 
the compressibility modulus K{ng), the effective nucleon 
mass m^(no) = nivr(l — /□), and the symmetry energy 
Esyai = ^{d'^E /The system of linear 
equations relating the saturation parameters and the pa¬ 
rameters of the model can be found, e.g., in Ref. 0. In 
our study we use as a reference the set of parameters: 
no = 0.16 fm“^ and Eg = —16 MeV and 


iF(no) = 250 MeV, E^y^ = 30 MeV, fg = 0.2. (5) 

The parameters of our NLW model are = 196.343, 
Cl = 90.7682, Cl = 88.7261, b = 8.94551 • and 
c = 7.70766 • 10“^, the nucleon mass is fixed as ttin = 
938 MeV. 

The composition of the NS matter being in the /3- 
equilibrium is governed by the relation between the lep¬ 
ton and nucleon chemical potentials Pe = Pp = Pn — Pp, 
where pn,p = dE/dun^p and by the charge neutrality 
condition requiring that the charge of protons is com¬ 
pensated for by the charges of electrons and muons, i.e., 
np = He + n^, with ni = 9{pf — mf){pf — TO/)^/^/( 37 r^) 
and mi standing for the mass of the lepton 1. Solution 
of this system of equations determines nucleon and lep¬ 
ton densities as functions of the total nucleon density 

n = Up + Un- 

The total energy density and pressure, E and P, for 
the NS matter are obtained after adding the kinetic con¬ 
tributions from leptons to the nucleon energy density and 
the pressure given by Eqs. ([3]) and (0]). With P — P{E) 
at hand, the NS masses and radii are calculated with the 
help of the Tolman-Oppenheimer-Volkoff (TOV) equa¬ 
tion. In the crust region we match smoothly our EoS 


with the Baym-Pethick-Sutherland (BPS) EoS [2^. The 
limiting mass of a NS for the present choice of parameters 
([S]) is equal to 1.92 Mq, which falls below the experimen¬ 
tal limit. 

Now we propose a simple method for controlling the 
dependence of the scalar field / on the density that, as we 
demonstrate later, allows us to stiffen the EoS at densities 
higher than a certain chosen density n, > ng. To be 
specific let us consider isospin symmetric matter (ISM). 
The function f(n) depends on the derivative of the scalar- 
field potential U'{f), see Eq. ([2|). Hence by changing 

U{f)^U{f) = U{f) + AU{f) (6) 

with an appropriately chosen AU{f) we may quench the 
fin) growth at densities n ^ n^. The requirement for 
the Uif) can be easily educed from the derivative of 
the scalar field over the nucleon density following from 

Eq. © 

(V ^_ 2idns/dn) _ 

dn m],Cy + U"if)/mN - 2(5ns/5/) ’ 

where U"if) stands for the second derivative of the po¬ 
tential U with respect to /. The partial derivatives of 
the scalar density are equal to 

PF 

dns dns f m^p'^dp/n'^ 

dn - 2 v'p^ + m;/’ df - J (p 2 + ^ ^ 

From Eq. © we see that although d//dn is always pos¬ 
itive [provided U''if) > 0 ] it can be made as small as 
required if one takes U"{f) sufficiently large for / in the 
vicinity of a value /core determined by the choice of the 
density n*. According to Eq. (© the requirement of a 
slow growth of the function /(n), ndf /dn <§; 1 translates 
into the condition 

AU”if) » mNm%, n/^jpl + m*^, (9) 

which is insensitive to the particular choice of the original 
potential Uif). 

Now let us give two examples of the appropriate de¬ 
pendence of AUif). We will refer to such an extension as 
the “u-cut scheme” and denote the model as the NLWcut 
model. The potential 

AU if) = a ln[l -h exp(/3(/ - /s.core))] (10) 

allows / to grow above the /s.core- This could be in¬ 
terpreted as an inclusion of the excluded volume effects 
with a soft core in the scalar baryon density. The max¬ 
imum of AU"if) in Eq. (ITUl) is realized for / = /s.core 
and is equal to a/3^/4. Hence condition © means 
/iAm\i) ^ 1.5 x 10“^(n/no). For numerical illus¬ 
trations we use the potential m with the parameters 
a = m\. = 4.822 x 10~'^ mff and /3 = 120 and vary the 
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FIG. 1: Left panel: the scalar field potential U{f) in the NLW 
model and U{f) in the NLWcut model for various values Ca, 
as functions of the scalar field parameter /. The symbols 
indicate the values of [/(n); U{f(n)), U{f{n))] for the density 
n of the ISM varying from no to lOno with the Ino step. Right 
panel: the nucleon effective mass in the ISM as a function of 
the nucleon density for the same models. 

value /s.core- Thus, we have /(Am%) ~ 1.7, and in¬ 
equality dH) is fulfilled. If we choose a singular potential 
diverging at / = /h.core, 

AC/(/)=a[<5//(/h,eore-/)]'^ (H) 

the function /(n) will never exceed the value /h.core- This 
could be interpreted as an inclusion of the excluded vol¬ 
ume effects with a hard core in the scalar baryon density. 
For densities n < lOno the same dependence /(n) as for 
the potential can be recovered with the potential in 
the form (fTTll with parameters a = 0.809 (3 = 11.138, 

and /h.core = /s.core + Sf with Sf = 0.190. 

Applying the method to models with different values 
of /o it is convenient to present 

/s.core = /o + C<r(l-/o). (12) 

For illustration we use the values = 0.2,0.3, 0.4, which 
correspond to /s.core ~ 0.36, 0.44, and 0.52 for /o = 0.2. 

The potentials U and U given by Eq. (Unj are shown 
on the left panel in Fig. [I]as a function of /. Symbols 
mark the values of [/(n); U{f{n))] and [/(n); U{f{n))] for 
various densities n. We see that indeed the symbols pile 
up in the region where the potential U has a maximal 
rise. On the right panel in Fig. [1] we show the effective 
nucleon mass m^(n) = miv(l — f(n)) as a function of 
the nucleon density, where /(n) follows as a solution of 
Eq. ([2]). We see that the function m^(n) flattens off for 
densities n > n* ~ 1.9 no, 3.0 no, and 4.0 no correspond¬ 
ing to Ccr = 0.2, 0.3, and 0.4, respectively. The closer the 
/s.core value is to /o, i.e. the smaller is, the smaller 
the density is, at which the NLWcut model starts to de¬ 
viate from the original NLW model. The NLW model 
is recovered from the NLWcut model for Co- ~ 1, when 
U{f) almost coincides with 17(/) for 0 < / 1 for all 

densities. We stress that all saturation properties of the 
NLWcut models remain the same as for the original NLW 



FIG. 2: Left panel: Total pressure P for the NLW and NLW¬ 
cut models as a function of the nucleon density in the ISM for 
various values Co-; the insertion demonstrates the binding en¬ 
ergy per nucleon. Two shadowed areas show the constraints 
on the pressure from the par ticle flow and the kaon production 
data in HIC extracted in 1^.1^. Bold dots show the extrapo¬ 
lation of the pressure consistent with GMR [^. Right panel: 
The NS mass and radius as functions of the central density for 
our models. The hatched band denotes the uncertainty in the 
value of the maximum measured NS mass of (2.01±0.04) Mq. 


model, since by construction AU{no) is chosen to be a 
tiny value. 

The influence of the quenching of the nucleon-mass de¬ 
crease on an EoS is illustrated in the left panel in Fig. [5J 
where we show the total pressure, Eq. o, and the bind¬ 
ing energy per nucleon as functions of the nucleon density 
in the ISM. We observe that the replacement (nni) leads 
to a sizable increase of the pressure for n > n* because of 
a reduction of the negative contribution to the pressure 
coming from the explicitly tr-field-dependent term, 
which overcomes a slight decrease of the effective kinetic 
term Pkin owing to a relative increase of the nucleon mass 
m*j^. Constraints on the pressure extracted from the par¬ 
ticle flow and the kaon production in heavy-ion collisions 
(HICs) [ 21 I, [23 are shown by the shaded areas and from 
analysis of the data on the giant monopole resonance 
(GMR) [i^l, by the bold dots. We see that the EoSs in 
the NLWcut model with the input parameters ([5]) satisfy 
these constraints for Co- > 0.2, exceeding the upper limit 
of the HIC constraints only slightly for Ca = 0.2 and in a 
narrow interval of densities, 2no <n < 2.6 uq. 

Different from the pressure, the energy-density and 
corresponding binding energy increase more weakly af¬ 
ter the inclusion of the new term 113 to the potential. 
This occurs because an increase in the effective kinetic 
energy Akin density is almost fully compensated for by a 
decrease in the explicitly cr-field dependent contribution 
Per. As a consequence, the net effect of the replacements 
(O and (fTUl) is a stiffening of the EoS P{E) for densities 
n > n,:. This result holds also for pure neutron matter 
and NS matter since the symmetry energy is not influ¬ 
enced by the replacement ®. 

The stiffening of the EoS is reflected in an increase 
in the NS masses as functions of the central density, 
which are shown in the right panel in Fig. [5] for our 
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FIG. 3: Contour plot of the maximum NS mass as a function 
of m%{no) and K{no) for the NLWcut models with Co- = 1, 
simulating the case of the original NLW model (left panel) 
and Cct = 0.3 (right panel). Contour lines are shown with the 
mass step 0.05 Mq. 


models together with the experimental constraint M = 
(2.01 ± 0.04) Mq. The proposed a-cut scheme allows us 
to shift the maximal NS mass from 1.92 Mq for the orig¬ 
inal NLW model to 1.96 Mq, 2.03 Mq, and 2.12 Mq for 
Co- = 0.4, 0.3, and 0.2, respectively. Thus, using our mod¬ 
ification one is able to fit the experimental constraint 
on the maximum NS mass, provided the value /s.core 
is chosen low enough (corresponding to Ca = 0.3 and 
0.2 in our example). Also, the NS radius R is shown 
in this panel. As one may see, R changes only slightly 
with Cct increasing with a decrease of Co-. For example, 
within the NLWcut models with Cg- > 0.4 a NS with 
the mass M = 1.4Mq has radius R = 12.7km, whereas 
R = 12.9 km for Ca = 0.2. 

It is instructive to study the dependence of the max¬ 
imum NS mass on the input values m^(no) and Kino) 
for different cut parameters Cn. In Fig. |3]we show the 
contour plot of the maximum NS mass for the origi¬ 
nal NLW model and the NLWcut model with Ca = 0.3. 
We see that the application of the cr-cut scheme removes 
completely the dependence of Mmax on the compressibil¬ 
ity modulus K{no) for m*^{no)/mN < 0.8, the contour 
lines on the right panel in Fig. [3] are almost parallel to 
the Ar('no) axis. The gain in M^ax proves to be larger, 
the smaller the ratio m)(^(no)/mjv is, reaching, e.g., at 
m%{no)/mN — 0.65, the value of 0.2 Mq for Ca = 0.3 
and 0.35 Mq for = 0.2. 

In this work we leave aside the question about hy- 
peronization of the NS matter at high densities, which 
would soften the EoS and decrease the maximum avail¬ 
able NS mass. A typical decrease of the maximum NS 
mass owing to hyperonization is rather large, typically 
~ 0.3-0.4Mq, so that the maximum calculated NS mass 
may become below the maximum measured NS mass (the 
so-called “hyperon puzzle” [1^). If included within the 
NLW model that we exploit here, the hyperons would ap¬ 
pear at the density 2.73 no, a value that changes weakly, if 
the (T-cut scheme is applied, decreasing with an decrease 
of Co- up to 2.56 no for = 0.2. A possible mechanism to 



FIG. 4: Same as Fig. [5] but for the FSUGold model. 


overcome the problem within the cr-cut scheme presented 
in this work is the use of a smaller values of m)(f(no) 
and/or Ca, see Fig. [3] However, in this case the problem 
is solved at the price of a possible violation of the HIC 
constraints on the pressure in a certain density interval. 
Other mechanisms to solve the puzzle might be an alter¬ 
native choice of the coupling constants, see [MLjfield- (or 
density)-dependent parameters of the model H, [l^, or 
something else. 

As a further illustration we apply the cr-cut scheme 
to the FSUGold model 0, which proved to be fine 
tuned to describe the properties of finite nuclei. This 
model includes the vector meson self-interactions, which 
are described by the interaction Lagrangian: £int = 
/4\ + with additional pa¬ 

rameters C and A„ determined as in Q. We modify the 
model by adding the scalar-field potential (nni with the 
same values of a and P and consider Ca = 0.2, 0.3, 0.4. 

The left panel in Fig. |T] shows the binding energy per 
nucleon and the pressure in the ISM for the FSUGold 
model and its cr-cut extension. The energy and pressure 
for the cr-cut model begin to differ significantly from the 
FSUGold model at densities n ~ 1.8 no for Ccr = 0.3 
and n ~ 1.4 no for Co- = 0.2. The saturation proper¬ 
ties of the FSUGold model are almost not affected for 
Ca £ 0.2. However, for Ca = 0.2 the pressure exceeds 
the HIG and GMR constraints at densities n < 3no. For 
Ca = 0.3 the constraints are satisfied except for a narrow 
density interval 1.7no < n < 2.7no. In the right panel 
in Fig. m we show the NS mass and radius as functions 
of the central density for the original FSUGold model 
and the FSUGold model with the scalar-field potential 
modified as in Eqs. m and for various values of Ca. 
The application of the cr-cut scheme allows us to increase 
the maximum NS mass from the original 1.72 Mq up to 
1.89 Mq, 1.98 Mq and 2.09 Mq for Ca = 0.4, 0.3, and 0.2, 
respectively. So, the model starts to satisfy the heavy 
NS mass constraint for Ca £ 0.3. Thus the re-tuning of 
the ESUGold model performed in (2^ would be unnec¬ 
essary, if one exploited the a-cut scheme. We note that 
the increase of the maximum NS mass is larger in the 
FSUGold model than in the NLW model with parame¬ 
ters this can be related to the choice of the smaller 
value of m^(no) = 0.61 niM in the former model. The NS 
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radius for the FSUGold model is more sensitive to the cr 
cut than for the NLW model, increasing for M = IAMq 
from 12.0 to 12.3, 12.7, and 13.3 km for Ca = 0.4, 0.3, 
and 0.2, respectively. 

In conclusion, we proposed a simple method of mak¬ 
ing the EoS obtained in a relativistic mean-field model 
stiffer at densities larger than some chosen value n*. The 
strategy is to quench the growth of the scalar field at 
the value /s.core- This can be achieved by adding to the 
energy-density a functional of the scalar field, which is 
vanishingly small for / ~ /o < /s.core and is rapidly 
increased for / ~ /s.core, so that its second derivative 
becomes sufficiently large. Then, the self-consistent so¬ 
lution of the corresponding equation of motion does not 
let the function /(n) to grow significantly above /s.core 
within a broad interval of densities, for n > n*. Herewith 
the EoS remains unaltered for densities n ~ no but stiff¬ 
ens significantly for n > n*, which results in an increase 
of the limiting value of the neutron star mass. We demon¬ 
strated the work of this method at hand of the standard 
non-linear Walecka model. Then we applied the method 
to the FSUGold EoS and showed that it is possible to 


increase the maximum neutron star mass bringing it in 
agreement with modern astrophysical constraints, with¬ 
out changing the EoS for densities below n ~ no. 

The value /s.core might be associated with a manifes¬ 
tation of the nuclear core in the medium and for / in 
the vicinity of /s.core one may expect occurrence of the 
deconfinement phase transition. 
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